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ABSTRACT 

Several new developments in the calculation and interpretation of hadron density- 
density correlation functions are presented. The asymptotic behavior of correlation func- 
tions is determined from a tree diagram path integral. A method is developed to use this 
behavior to correct for leading image contributions on a finite periodic spatial lattice and 
to correct for the finite temporal extent of the lattice. Equal time correlation functions are 
shown to determine a sum of the ground state rms radius plus a polarization correction, 
and it is shown how to extract the hadron polarizability from unequal time correlation 
functions. Image-corrected correlation functions calculated in quenched lattice QCD are 
presented and used to determine the size of the pion and nucleon. 



-1- 



I. INTRODUCTION 

One of the primary motivations for solving QCD on a lattice is to understand the 
structure of hadrons. Hence, in addition to reproducing the hadron spectrum and ex- 
perimentally measured form factors, it is important to explore the spatial distributions 
of quarks and the correlations between them in the ground states of hadrons. As in 
other strongly interacting many-body systems, two-body correlation functions are a natu- 
ral starting point, with the simplest being the density-density correlation function 



where \h) denotes a hadronic ground state and u and d represent field operators for up 
and down quarks respectively. Although these correlation functions have been studied 
in a number of works [1-4], several important practical and conceptual problems were 
previously unresolved and errors were made in interpreting correlation functions. Hence, 
this work presents a number of new results relevant to the calculation and interpretation 
of density-density correlation functions. 

One practical problem surrounding all lattice QCD calculations is the irreconcilable 
conflict with finite computer resources between reducing the lattice spacing (and thus 
accurately approximating the continuum limit) and increasing the physical volume of the 
lattice (and thus eliminating finite volume corrections). Nowhere is the confiict more severe 
than in the calculation of density-density correlation functions in hadrons. Conventional 
hadron wave functions, either gauge fixed in Coulomb or Landau gauge or made gauge 
invariant by a line of flux e'f'^^ fall off much more rapidly tharr quark correlation 
functions because they contain additional suppression factors at large separation arising 
from the small overlap between the gluon wave functional in the hadron and in the vacuum 
[5] . Typically the spatial extent of these wave functions is half that of the density-density 
correlation functions, so that lattices in common use which are adequate for wave functions 
are much too small for correlation functions [3] . 




(1.1) 
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In section 2, we therefore take a fresh look at the general problem of relating corre- 
lation functions in a periodically replicated array of hadrons, corresponding to periodic 
boundary conditions on a spatial lattice, to those in an isolated hadron. We begin with a 
pedagogical example in one spatial dimension and show how to correct for the effects of 
periodic images and generalize to the physical case of three space dimensions. A crucial in- 
gredient in obtaining the properties of an isolated hadron is understanding the asymptotic 
decay of the correlation functions. Since the behavior is more complicated than simple 
exponential decay governed by the rho mass, in section 3 we construct a tree diagram path 
integral which incorporates the essential physics. This path integral accurately describes 
the asymptotic behavior and provides a means to remove finite spatial and temporal volume 
effects. It also explains an apparent paradox wherein two seemingly equivalent methods 
of calculating correlation functions yielded different asymptotic behavior [6] . 

Given the possibility of accurately determining the density-density correlation function 
for an isolated hadron, it is important to understand clearly its physical content and to 
explore the full range of physical information which can be extracted from it. For this 
reason, in section IV we reexamine the physical content of the density-density correlation 
functions at equal and unequal time in the nonrelativistic limit, which is particularly simple 
to understand. We show that the rms radius of the equal time correlation function is not 
the rms radius of the hadron, as previously claimed, but rather is the sum of the rms radius 
and a dipole-dipolc polarization term. Furthermore, we show how to extract the hadron 
polarizability from integrating the correlation function over relative time, and generalize 
the results to the relativistic case and more general currents. 

Finally, having developed the general formalism and methodology, the results of 
quenched lattice calculations for density-density correlation functions for the n, p and nu- 
cleon are presented and discussed in section V. The conclusions and outlook are discussed 
in section VI. 
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II. PERIODIC IMAGES IN LATTICE GAUGE THEORIES 

Intuitively, it is clear that when we calculate hadron correlation functions in a fi- 
nite spatial volume with periodic boundary conditions, we are not considering an isolated 
hadron but rather a periodically replicated array of hadrons. Consequently, periodicity ef- 
fects and interactions between neighboring hadrons will occur and we either must discover 
how to correct for these effects or calculate with a volume so large that these effects are 
negligible. 

The primary result we will establish in this section is that the density-density cor- 
relation function for a gauge theory on a periodic lattice of spatial dimension L has the 
form 



where poif) is an accurate approximation to the free hadron result p(r) for r < L and 
differs due to interactions with periodic images only for r^L. This result embodies two 
physically distinct effects. The dominant effect of periodicity is just the summation of the 
tails of po from neighboring images, and this effect amounts to a factor of 2 correction 



and p arising from interactions with periodic images, and it first appears at distance L. 
Fortunately, since the dominant effect is the sum of tails from images, it can be removed 
by subtracting these image contributions, leaving only the residual discrepancy between 
Po and p beyond r ~ L which will be negligible in cases of practical interest. 

To understand how the form (2.1) arises in a gauge theory and how po may be obtained 
by subtracting image contributions, it is useful to begin by considering the simple example 
of QED coupled to non-relativistic particles in 1 + 1 dimensions on a circle. 

Positronium on a circle 

In the continuum, using the gauge d^A^ = 0, the Hamiltonian for a non-relativistic 
e'^e~ pair in 1 + 1 dimensions on a circle with circumference L is 




(2.1) 
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at distance along the Cartesian axes. The second effect is the discrepancy between po 
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where V{x) is the periodic potential 
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The Schrodinger equation corresponding to H, Eq. (2.2), must be solved with appropriate 
boundary conditions on the wavefunction. Here we will use periodic boundary conditions 
in a;i and X2 

Ip {Xi + L, X2, Aa;) ='lp{xi,X2 + L, A^) = ip {xi,X2, A^) (2.4) 

although other boundary conditions are also possible. The center-of-mass motion separates 
in Eq. (2.2), i.e. 

2L dAl Am m ^ ^ ^ ^ 

where P = pi + P2, P = (pi — P2)/2, x = xi — X2- Notice that Eq. (2.4) implies period- 
icity in the center-of-mass and the relative coordinate separately. The eigenstates of H 
corresponding to P = can be expressed in the form 

M^, Ax) = Yl e'^"""^^''^"V^fc(a; - nL) . (2.6) 

n 

where is a solution of the equation 

EkMx) = + ^i^)) Mx) , (2.7) 



with 



p2 

V{x) = V{x) + —x'^ . (2.8) 



Note that i(^{x, Ax) is defined only for —L < x < L (with periodic or quasiperiodic bound- 
ary conditions), whereas V'(x) is defined for — oo < x < oo (with vanishing boundary 
conditions). This difference in domain allows a function of one variable to carry the same 
information as a function of two variables. 



For X = nL + Ax, where n e IN and < Aa; < L, the effective potential V can be 
rewritten in the form 

V{x) = ^{nLf + {2n + l)^Ax , (2.9) 

which allows a simple physical interpretation. Consider the case in which the e'^ moves 
around the circle n times while the e~ is fixed. The electric flux from the e"*" to the 
e~ ends up in a configuration where it is wrapped around the circle n times, since there 
are no terms in the Hamiltonian which would allow the gauge field to make a transition 
to the energetically lowest unwrapped configuration. The electric field energy density 
is proportional to the electric field strength squared — the latter being (n + 1) • e for 
< X < Ax and ne for Ax < x < L in the above example. Hence 

^gauge = l{^^- [(^ + + iL- Ax)[ne]^'j (2.10) 

which agrees with Eq. (2.9). 

Note that for —L < x < L the effective potential (which determines the excitation 
spectrum of H) agrees with the potential on a line, i.e. V{x) = ^\x\. In contrast, if we 
had omitted the degree of freedom, we would have obtained {\x\ < L) 

V{x) = ^\x\ - , (periodic) (2.11) 

2 

which is a much worse approximation to \\x\ than Eq. (2.9). Although the Ay; degree of 
freedom "freezes out" for L ^ oo (see Eq. (2.5)), its presence improves the approach to 
the infinite volume result for finite L. 

For I a; I > L, the interaction is modified by the periodic boundary condition. In a sense 
the positronium interacts with its own periodic replicas. However, in most cases, this has 
only an exponentially suppressed effect on energy eigenvalues and the wavefunctions, since 
the wave functions fall off exponentially for large x. 

For ground state correlations another finite size effect is of much more importance. 
Consider 

pj,{x) = Jdy {k\5{xi - y)S{x2 - x - y)\k) , (2.12) 
-6- 



where \k) denotes the k^^ positronium state. Using Eq. (2.6), 
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where 



pk{x) = ^k{x) 
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(2.14) 



The essential result is that the interference terms (oc ipni^m^ n ^ m) are removed by the 
integration over A.* 

Equation (2.13) shows explicitly how the general result, Eq. (2.1), arises in 1 + 1 
dimensions, gives a precise definition of po as p/s, and displays the role the gauge field 
plays in obtaining this result. The longitudinal gauge field A^, which is the only quantum 
mechanical degree of freedom associated with the gauge field in QEDi+i, plays two essential 
roles. First, it removes all cross terms so only the square of the wave function |7^fc(a;)p 
appears in the final result. Secondly, it replaces the periodic potential, Eq. (2.11), by the 
linear potential in the range —L < x < L, so that i/jki^) satisfies the correct Schrodinger 
equation in this entire range. The role of Ax has also been emphasized by a number 
of authors in the context of the Schwinger model [7], where it turns out to be crucial in 
providing a proper description of the anomaly in the divergence of the axial vector current. 

A numerical example in 1+1 dimensions 

To show how the density for a free bound state, po{x), may be extracted from the pe- 
riodic sum Pperiodic(a;) in Eq. (2.1), we now consider a numerical example. For convenience, 

* To avoid irrelevant details, we have not been explicit about normalization and the 
5(0)-term (m = n in Eq. (2.13)) arising because we work in the continuum. In a lattice 
calculation, the gauge field integration is compact. 
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and to establish contact with the lattice formulation of QCD3+1, we use a discrete lowest 
order difference approximation to the Hamiltonian, Eq. (2.7). Note that beginning with 
QEDi_|_i lattice gauge theory and gauge fixing analogously yields an equivalent form with 
one quantum mechanical gauge degree of freedom remaining [8,9]. The general conclusion 
will be independent of the discretization. 

To emphasize the image corrections, we have chosen the example shown in Fig. 1, in 
which there is greater overlap between the bound state and its periodic images than in a 
typical lattice QCD calculation. Note that the periodic solution , denoted by the heavy 
solid curve, falls to less than half of its peak value. This case has N = 10 lattice sites 
separated by spacing a, with mass ma = 0.3 and charge ea — 0.3. With these parameters, 
the exact solution on the open line, shown by the light solid curve, has size — 3.6a. 

Our general strategy for image corrections is to use the known asymptotic decay of the 
infinite domain solution both to subtract the tails of the first images in the fundamental 
unit cell and to approximate the solution outside the unit cell. Thus, in any dimension, 
assuming the asymptotic form Pas(r) known, the image corrected density Pcor(r) 
defined 

Pcor(r) = Pas(^) (2.15a) 

for r outside the fundamental unit cell and 

Pcor (r ) = Pperiodic (r ) - ^ Pas (r + nL ) (2 . 1 5b) 

for f inside the first unit cell. 

In our one-dimensional example, we make use of the known asymptotic behavior for 
Airy functions, so that as a — > 0, the density has the form 

Pas(a^) ~ exp (^-Ax^/^j , where X=-J—- . (2.16) 

In the continuum limit, with the above parameters, A = 0.155. For finite lattice spacing, A 
is somewhat smaller. The normalization of Pas(a^) is determined by the symmetry condition 

Pas(^) = ^ Pperiodic(^) (2.17) 

-8- 



Although, by construction, Pcor(a^) and its first derivative are continuous at a; = its 
second derivative is not, both because of finite a corrections to A and because p{x) is not yet 
asymptotic at a; = Especially when one is also interested in the Fourier transform, 
it is desirable to parameterize the image corrections such that the second derivative is 
continuous, and there are several possibilities. One possibility is simply to define A such 
that Pcor(^c) also has a continuous second derivative. Alternatively, one may make the 
ansatz 

p^{x) ~ (a + hx) exp (^-Aa;^/^) (2.18) 

and determine h/a by the same criterion. Since there are polynomial corrections to 
Eq. (2.16), this form is clearly reasonable. The results are essentially the same, and we 
only show calculations using the ansatz (2.18). 

The result of correcting the periodic solution using Eq. (2.15) with the asymptotic 
form Eq. (2.18) with A = 0.13 is shown by the dashed line in Fig. 1. This corrected result 
agrees quite well with the exact solution on the infinite domain, shown by the light solid 
line, especially considering the high degree of overlap in our deliberately chosen worst-case 
example. Note that the reconstruction of the density is not extremely sensitive to the 
precise value of A. Although the discrete solution on the open line yields A w 0.145, the 
value A = 0.13 we have chosen and A = 0.15 used in calculations which are not shown here 
yield very similar results. 

The Fourier transforms of the densities in Fig. 1 are shown in Fig. 2. We will show 
in Section IV that they are the sums of squares of elastic and inelastic form factors. The 
Fourier transform of the periodic density is only defined at the discrete lattice momenta 
qk = kiT / L and is indicated by the solid dots in Fig. 2. An equivalent way of rephrasing 
our image correction problem is to ask whether we can use our knowledge of the form of 
the density, Eq. (2.1) and the asymptotic behavior to reconstruct the Fourier transform 
for all qa. Clearly, the curvature at the origin, for example, is related to the rms radius 
of the correlation function and must depend sensitively on the details of the surface in a 
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way which cannot be guessed by naive polynomial interpolation of the solid points. One 
observes that the Fourier transform of the corrected density from Fig. 1, denoted by the 
solid curve in part (a), agrees well with the Fourier transform of the exact result on the 
open line. We should emphasize that the level of agreement shown in part (a) of Fig. 2 
and in Fig. 1 is not unique to the particular set of parameters we have chosen, and that 
similar results were obtained for a wide range of parameters N, ma, and ea. 

The importance of making the second derivative of the corrected density continuous 
is shown by comparing curves (a) and (b) of Fig. 2. The solid curve (b) shows the result of 
using Eq. (2.16) for the asymptotic tail with the same value A = 0.13 but not introducing 
an additional linear term as in Eq. (2.18) which can be used to make the second derivative 
continuous. The substantial oscillations around the exact result show how large the effects 
from the discontinuity become at intermediate values of qa. One should note that Eq. (2.1) 
has the property that at the discrete lattice momenta qk = kir/L, the Fourier transforms of 
Pperiodic and po circ identical, so that to the extent that second and higher image corrections 
are negligible, both the solid and dashed lines must go through the solid dots. 

In practical lattice calculations, the periodic signal pk{x) in Eq. (2.13) is not a smooth 
function but rather is contaminated with statistical noise. This adds another level of 
complication to the reconstruction of Pcot{x). To study this case, in part (c) of Fig. 2 we 
have considered an ensemble of 10 data sets obtained by adding randomly distributed errors 
with mean fractional deviation 1% to the periodic density at each lattice site. As before, 
for the matching conditions we use continuity of p and the second derivative. Instead 
of using a three-point formula, both are now computed using a five-point formula which 
automatically provides some averaging over the statistical fluctuations. One observes that 
for q < tt/L the resulting form factor turns out to be quite stable, and one has as before 
an accurate determination of the rms radius. However, as expected, for larger values of 
q at which the Fourier transform has fallen by two orders of magnitude, the 1% errors 
introduce uncertainties comparable to the signal. In this region, the uncertainty is well 
represented by the variance of the solid dots coming from the periodic solution and no 
further information can be obtained by image corrections. 



Image Corrections in QCDg,,.]^ 

Fortunately, although one cannot carry through an explicit solution as in the case 
of QEDj^^.]^, Eq. (2.1) also applies to QCDg^.^ with periodic boundary conditions in the 
spatial direction. The effects of images may be understood by considering all possible 
contractions of the field operators occurring in all periodic replicas of the density operators 
and all periodic replicas of the sources. 

Representative contractions relevant to meson density-density correlation functions 
are sketched in Fig. 3, where one periodic spatial dimension x is shown on a cylinder and 
the other two spatial dimensions are suppressed. The circumference of the cylinder is L, and 
the separation between the two density operators on the midplane of the cylinder measured 
on the front of the cylinder is denoted r and around the back of the cylinder is \r — L\. The 
upper sketch shows the case in which all contractions occur within the fundamental unit 
cell and yields the physical result in the infinite volume limit. The middle sketch shows 
a typical contraction in which the meson is created in the fundamental unit cell on the 
left creating the wave function ipt^ir) at the midplane whereas on the right, a propagator 
from the first periodic image creates the wave function 'il)Q{r — L) at the mid plane. For 
a periodic wave function which could be written in the form '0(r) = V'o(^+ ^L), this 
contraction would thus correspond to a cross term of the form '^^{r -\- nL) il;o(f + mL). 
Clearly, if this term were non- vanishing, one would never obtain the form Eq. (2.1) which 
is just the sum of diagonal terms. However, since the total propagator in the middle 
sketch has the topology of a spatial Polyakov line, it therefore vanishes in the confining 
phase for a sufficiently large lattice. This is the physical reason no cross terms occur in 
gauge theories in any dimension. Finally, the bottom sketch shows a typical contraction 
in which the meson is created with propagators from the first periodic image on the left 
and on the right, corresponding to the diagonal term i/^oii^ ~ L) i(^o{r — L), giving rise to a 
diagonal contribution to the density of the form p(r + nL). Extending this argument to 
all contractions in 3+1 dimensions, one can see that the Polyakov line argument removes 
all cross terms, and that the result must have the form of Eq. (2.1). 



Thus, in principle, the procedure for correcting images in higher dimensions is com- 
pletely analogous to that in two dimensions. However, from the practical point of view 
there is an additional complication since the boundary of a cube in D > 1 dimensions is 
a D — 1-dimensional extended object and it would be very difficult to impose continuity 
everywhere on the boundary. Therefore, it is useful to consider more general matching 
conditions to replace the continuity requirement. 

To this end, it is important to note that there are only a few critical points on the 
boundary where it is essential to achieve continuity. This can be seen by calculating the 
angular averaged form factor 

p{q)= [ dVLq~p{q)o, H drr^-'^^^p{r) , (2.19) 
J Jo 

where 

p(r) = J dflrPir) . (2.20) 

Obviously, if one wants to avoid unphysical oscillations in p{q), one should keep p(r) as 
smooth as possible. In order to see what conditions this implies for p(r), consider the 
extreme case of a step function in two dimensions {r — {x,y)) 

fO \x\>aoT\y\>a 
' \l |a;| < a and |y| < a ^ ^ 

yielding, after taking the angular average, 

1 < r < a 

4 r 

1 sin~^ — a<r<a-^/2 (2.22) 

TT 2a ^ ' 

r > a- \/2 

Note that although p(f) is discontinuous everywhere on the boundary of the square, p(r) 
and its derivatives are continuous everywhere except at two critical points: r = a and 
r = a/2 a, where p'{r) is discontinuous. Physically, this corresponds to that fact that 
when one increases r and the sphere |r| = r intersects with the corner or the center of 
the faces of the square (Eq. (2.22)) then there is a sudden change in the fraction of the 




sphere contained inside the square, reflected in a discontinuous derivative. Turning now to 
three dimensions, we see that if we need to match solutions inside and outside the cubic 
unit ceU where discontinuities can occur only at the boundary, discontinuities in p(r) or its 
derivative will arise only if p(r) is not continuous in the corners, the center of the edges, 
or the center of the faces of the first unit cell. Hence, unphysical oscillations in the form 
factor at large q are avoided most efficiently if the parameters of the asymptotic density, 
such as in Eq. (2.18), are determined such that the discontinuities around these critical 
points are minimized. 

In the practical case of lattice calculations in which the density measured at every 
point has statistical errors, a different strategy is required, and an iterative self-consistent 
method to determine the optimal density is described in section V. 

III. THE ASYMPTOTIC BEHAVIOR OF DENSITY-DENSITY CORRE- 
LATIONS 

In the last section we showed how one can use the knowledge of the asymptotic 
behavior of the density for an isolated hadron, po(^), to correct for images in the density 
of a periodic system, related to po by 

Pperiodic(r) = ^ Po (r + nL) . (3.1) 

n 

Hence, because our ultimate objective is to study lattice measurements of physical 
correlation functions, we will show in this section how to determine the behavior at large 
separation \y\ of the pion density correlation function: 

p{y, tiM) = 5^(7r|p"(:?, h)p''{x + y, ^2) \n) (3.2) 

X 

= Yl J](^^n| J+(rout, T) p"(f , ti)p''{x + y, t2) J{nn, 0)\nn) 

X n 

where the up quark density operator is p'^{x, t) =: u^^t 7° 'iJ'x,t '■ , p^ is the corresponding 
operator for down quarks, J{x, t) is the pion source ug^t 75 ug^t ■> {^n} denotes states 



with the quantum numbers of the vacuum, fin and fout denote the positions of the pion 
sources, and we work in Euchdean time. The analysis of the hadron correlation function 
is analogous. 

We will use the fact that at large distances, QCD is dominated by light hadrons to 
motivate a tree graph approximation to the asymptotic correlation function. Consider a 
typical time history for the pion correlation function p{y) shown in the left portion of 
Figure 4, where for simplicity we consider the equal time case ti = t2 = t. In the quenched 
approximation, one has a ud quark-antiquark pair connecting the sources J, J'^ and the 
density operators and p'^ as shown by wiggly lines, interacting via the exchange of gluons 
which are suppressed in the figure. Physically, at large Euclidian separation we expect the 
interacting quark-antiquark pairs to form meson ground states in the appropriate channels, 
so that the quark-level time history in the left sketch is replaced by the corresponding 
meson tree graph shown at the right. As long as each meson propagator is of sufficient 
length, this approximation should be accurate. Formally, the same result may be obtained 
by inserting complete sets of meson states between each of the operators and identifying 
the leading singularities [10]. In either case, since the sources couple to the n and the 
external currents couple to the p, the tree-graph approximation to the correlation function 



where for simplicity we use lattice propagators for scalar mesons [11] for the D's. In 
practice, to evaluate this expression on a lattice for comparison with lattice QCD results, 
it is preferable to Fourier transform to momentum space and perform the sums over discrete 
lattice momenta. 

A further approximation to the sum over tree graphs in Eq. (3.3) which provides insight 
into the asymmetric behavior is given by the stationary or classical approximation in which 
instead of summing over all joint positions ri and r2, we select only those joints in which 



Eq. (3.2) yields 




X ri,r2 



(3.3) 



X gpnn Dp {x - fi) gpn^Dp (x + y - ra) + ("1 ^ 2") 
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the product of propagators takes on maximal values. Replacing the lattice propagator by 
a simple continuum exponential D(r) — > e""^*", we obtain 



Pstationary(^) = J^iPp^nf Hiax e"^-!^^-^'"! e""^-!^^-^^! e"^^!^-^^! e-^'-l^+^-^^l e-^^'l^-*- 

(3.4) 

Here | | denotes the Euclidean distance and r-m = (0, 0), rout = (T, 0), and x = {ti,x)). 
The physical picture that results depends on the numerical value of the ratio m.j^/mp. 

The asymptotic behavior of pion correlation functions calculated in lattice QCD and 
in the tree-level approximation are compared in Fig. 5. Although the lattice calculations 
are discussed in detail later in Section 5, at this point the only information that is needed 
is the fact that calculations were performed at three values of the hopping parameter k 
corresponding to the hadron masses given in Table I. For the heaviest quark mass case, 
K2, TOtt is only 27% below rup and for the lightest case, K5, rriT^ is slightly above \mp. 

One important result seen in Figure 5 is the fact that beyond 1 fm, the tree-level 
approximation accurately describes the behavior of the lattice QCD results denoted by the 
solid curves. The tree-level results at selected points before image corrections are shown by 
the squares: Three interior points for which image corrections are negligible, are denoted by 
solid squares, and results at all other points which are subject to image effects are denoted 
by open squares. By symmetry, image contributions at the center of the faces, center of 
the edges, and corners of the first unit cell may be subtracted by dividing by 2, 4, and 8 
respectively, and the results so corrected are denoted by the solid octagons. Thus, all the 
solid symbols are free of image effects and are seen to agree extremely well with the solid 
lattice QCD curves as claimed. For the purpose of this comparison, the image- corrected 
tree level results were normalized to agree with the lattice results at r/a = 14. 

Another significant feature observed in Fig. 5 is the fact that in the region of 1-2 fm, 
the slope does not approach a fixed value given by a meson mass. Rather, it gradually 
changes in this region, taking on values intermediate between exponential decay with rUp 
and TTT-TT indicated by the dotted lines in the figure. At first, this may seem surprising. 



since because the current couples to the p, one might naively expect the asymptotic decay 
of the correlation function at large distances to decay as e~'^p'^ . 

The behavior of the asymptotic slope may be understood qualitatively using the sta- 
tionary approximation to the tree-diagram sum. A useful way to state the stationary 
condition for the joints ri and f2 in Fig. 4 is to connect all the vertices by classical strings, 
weight the length of each string segment by the corresponding meson mass, and minimize 
the total energy of these classical strings. In this approximation, it is simple to determine 
the asymptotic behavior in various regimes of interest. 

First, consider the case ttt-tt < |"T,p, which corresponds to the physical case and is 
also the simplest technically for the above model. The minimizing configuration is always 
obtained by shrinking the p-meson line to zero and connecting sources and currents with 
pions, that is, the joints coincide with the location of the currents. Thus, for r <K T, 
p{r) oc e""^'^'^ while for r » T, p(r) oc e"^"^'^'^ since the pions have to go back and forth. 
The case \rap < m.j^ < nip is similar except that it is energetically more favorable for 
the two pions in the t-channel to combine into a p-meson when the angle between them 
becomes smaller than the critical angle 

Combining both cases, we conclude that one should observe a transition from p oc e"'^''^ 
to p oc e~'^"^ where uih — min{2m7r, rup} as r increases from r <K T to r ::>> T. 

Physically, one should note that the transition from p oc e~'^'^'^ to p oc e""^"^ when 
r ~ T is a lattice artifact and that the true physical result corresponds to the limit T ^ oo. 
This limit is straightforward to evaluate using the tree level path integral result and thus 
provides an extremely useful extrapolation tool. In practice, one first calculates the tree 
level path integral for the actual lattice geometry as done in Fig. 5 to verify consistency 
with the lattice result. One then recalculates the tree level path integral for T — > oo, and 
uses the change in the asymptotic behavior as a correction to the finite lattice result. The 



effect of this correction is shown by the dashed hnes in Fig. 5 at K2 and K5, which show 
the correct physical asymptotic behavior in the hmit T ^ 00. 

In principle, the tree diagram path integral analysis is analogous for the p and nucleon. 
However, an interesting new feature arises in the case of the p. Naively, one would expect 
the relevant diagram to be of the form shown in Fig. 4 with all three tt's replaced by p's. 
However, the operator u^^u ~ — p^ does not project onto a physical p but rather 
creates a linear combination of a p and uj. The uj component may then couple to a p and 
TT, since the effective hadron theory contains a 7ru>p vertex of the form 

^Tvujp = gTvujpda^^fidj'Kpse"''^'^^ (3.6) 

Hence, there is also a diagram in which p's propagate from the sources to fi and r2, <^'s 
propagate from p" and p*^ to ri and r2, and a tt propagates between ri and r2- Although 
^TTLjp niay be suppressed by the derivative couplings, at sufficiently large |ri — r2|, the pion 
mass in the propagator will ultimately dominate the asymptotic decay, and we will observe 
this behavior in the lattice results presented in Section V. 

The tree level diagram model also allows one to understand the previously puzzling 
discrepancy between two calculations of density-density correlation functions in Ref. [6]. 
One calculation, which we will refer to as projected, integrated p{x, ^) p{x + y, f-) over 
all X to project onto zero momentum as in Eq. (3.2). The unprojected calculation used 
p(0, ^) p{x, -^) and relied upon the fact that non-zero momentum modes would automati- 
cally be suppressed by evolution for sufficiently large imaginary time T. 

The argument is simplest in the stationary-phase approximation to the tree level path 
integral, which we have verified is qualitatively similar to the full numerical integral. The 
dominant term in the sum Eq. (3.2) arises from x = — y/2. If we compare the resulting 
string configuration sketched in Fig. 6 (which maximizes Eq. (3.4)) with the one for x = 0, 
one observes an important difference. For small y, the minimal configuration arises in 
both cases from pure 7r-exchange. Since ttt-tt > mp/2, p-mesons start to develop in both 



cases for large y. However, in the unprojected case the relevant angles are smaller and the 
p-meson strings therefore develop earlier as shown in the sketch. Hence, we expect that 
the unprojected correlation function will approach the p slope earlier than the projected 
one, and this expectation is verified in the calculations shown in Fig. 7. This result agrees 
qualitatively with that presented in Ref. [6] and, we believe, explains the reason for the 
observed behavior. However, physically, we know that the approach to the p slope rather 
than the tt slope is an artifact due to the fact that the extent in T is too small. Thus, 
the most physical calculation is to perform the projection summation, which significantly 
suppresses the p contributions and then, in addition, correct for the effect of finite T by 
calculating the difference between the tree diagram path integral with finite T and T ^ oo 
as discussed in connection with the dashed curves in Fig. 5. 

IV. DETERMINATION OF HADRON PROPERTIES 

FROM CORRELATION FUNCTIONS ON A LATTICE 

In this section we will study density-density correlations on a lattice [12,2] 

P{y,tl,t2)'^= ^{hs \{u3^t2l^U3^t2) {dx+y,til^d$+y,ti)\ hs) , (4.1) 

X 

where \hs) is some superposition of states created by our lattice sources at to = and 
annihilated at = T. Inserting complete sets of states one finds 

g-(T-t2)p°^ 



Piy,h,t,)= J2 J2''''^*pM^o (4-2) 

hi,h2,n p,q '^a 

(/i2,p |u7°u| n,p + q) ^ {n,p + q {d'y^d] hi,p) ^ C^,^, 

Pn Phi 



X 



The kinematic factors = \Jp^ + M^, p° = ■\J{p + qf + have been introduced to 
insure that the states satisfy covariant normalization conditions in the continuum limit. 
Equation(4.2) is exact, but as it stands, is not very useful because it is too complicated. 

In addition to summing over all excitations n in the intermediate state, we have also 
summed over all excitations of the initial {hi) and final (/j.2) states as well as over all 



momenta to which the lattice sources couple, with amplitudes Cf. In this work we will 
use bag model sources [13] to produce a maximal overlap with the ground state hadrons. 
However, it is technically impractical to project such localized sources to a fixed momentum 
such as p = 0. Since localized sources imply an uncertainty in momentum space, one has 
to pay a price for physical sources by accepting a convolution over momenta. 

In many cases, Eq. (4.2) can be simplified considerably. Although the bag sources 
imply a convolution in momentum space, the large overlap with ground state hadrons 
allows one to drop the sum over initial and final state hadrons (provided the currents are 
not too close to the walls) and one is left with the correlation function in the ground state 
\h) 

X (h^p \u'y^u\ n,p+ q) ^ (n,p + q d\ h,p) e*^^ , 

Pn 



where Po = P° = yP^ + M^, t = t2 — ti. As discussed above, we are a priori not allowed 
to drop the summation over the momentum of the initial/final state. 

Let us first consider the simpler equal time case, t = 0. In this case, the second 
moment of the density-density-correlation is related to the size of the hadron. For example, 
for non-relativistic states one finds 

E^V (y, 0) 



(4.4) 



y 

where we have introduced the non-relativistic form factors 



5=0 



F}:^{q)^{h,p\u^'u\n,p + q)^{ 6 ' ^ , (4.5) 



q-df^^ + O (g2) n^h 



and similarly for the down quarks. 
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Inserting the low ^ expansion of the form factors (4.5) into (4.4) and separating the 
result into obvious ground state and polarization terms, one thus obtains 



pol 



(4.6) 



The origin of the dipole-dipole term, which has been omitted previously [2,4], becomes 
particularly clear if we rewrite the correlation function in the form 

2 
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h 



h 



1 



q- \ ru-rd 



q ■ r. 



n 



q ■ Td 



h )+0 



n ) {n\h) 



h 



it) 



(4.7) 



n } (n 



Td ■ q 
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where we have inserted complete sets of states in order to have only single-particle operators 
appearing. Due to the orthogonality of the states, only n = h contributes in the first two 
terms and these terms give rise to the ground state rms radius in Eq. (4.6). The last term, 
however, includes the dipole-dipole transition contributions to the correlation. 

Let us now estimate the importance of the d ■ d term for nonrelativistic quarks. From 



the positivity of ( (r„ ± r^)^ ) one finds the rather crude bounds 



pol 



gs 



(4.8) 



For more realistic estimates one has to make some model assumptions. For example, in a 
non-relativistic two-body problem with equal masses for up and down quarks one finds in 
the center-of-mass frame f„ = — and thus 



{K-7dr} = 2{Rl + Rl) , 



(4.9) 



so that, in terms of our definitions, 

('•'>po, = ■ 

More generally, defining Y,i<j{{i^i - fj)^) = N{N - 1) [(r^)gs + (r^)poi], one finds for an 
A/"-body bound state with equal masses and a symmetric coordinate space wavefunction 

In this non-relativistic estimate, for the nucleon (r^)^^^ = 1/2 (r^)^^ so that neglect of 
the polarization correction would lead to an overestimate of (r^)gs by 50% — certainly 
not a negligible correction. The decrease of the relative importance of {r'^)^^^ with the 
number of constituents, Eq. (4.11), can be understood easily. For N ^ oo the motion of 
two constituents becomes more and more uncorrelated and more accurately described by 
a mean field, so that (r^ • Vj) ^-^^ for z 7^ j. Thus, ^(r^ — r^)^^ (rf ) + (r|) for i 7^ j. 

The physical case involves significant corrections to the previous non-relativistic ar- 
gument. A sizable fraction of the mass can come from the gluons and the center-of-mass 
of the quarks need not coincide with the center-of-mass of the hadron. One can even 
imagine an extreme situation, in which the quarks are tightly bound together and move 
(together) around a large cloud of gluons — thereby yielding a large value for {i^'^)^^ while 
^(fj — Tj)^^ remains small — which provides a scenario where one approaches the lower 
bound in Eq. (4.8). In addition, relativistically, the notion of a center-of-mass is no longer 
appropriate, since the center-of-mass fluctuates. Although we are unable to provide a 
realistic estimate for QCD, the importance of the d- d-term in Eq. (4.6) should be evident. 

In a lattice calculation, one would prefer not to rely on model assumptions, so it is 
helpful to consider density-density correlation functions at unequal times [14] . Neglecting 
the motion of the center-of-mass for a moment, one finds in the non-relativistic case 

m ^ ^^^^ = «> + (rl) -'^Y.^-^^^-^^^'dlXn ■ (4.12) 



Thus the d ■ d-terms are exponentially suppressed for large t, which allows one not only 
to extract {r'^)+{r'^) in a model independent way but also to extract the oflF-diagonal 
elements of the polarizability tensor aud (where the rows and columns are labeled by flavor 
indices) from the approach to the asymptotic value by means of 

^ ->u ->d 

/ dt[R\t)-R\^)]=-2Y,p^='^o.ud . (4.13) 

A direct lattice measurement of hadron polarizabilities would be extremely interesting in 
view of efforts to measure them experimentally [15]. 

However, as we indicated already above, there are extra complications arising from the 
motion of the hadron between the two measurements as well as due to relativistic effects. 
With* 



'f + Ml + ^{p+qf + Ml + ^q^^ + O {f) n = h 



(4.14) 



q + O (q^) 



n ^ h 



where 



q'={^Jf + Ml-^(S, + qf + Ml^ - f = -f + + O (f) , (4.15) 

one finds for the second moment of the correlation function, after averaging over the angular 
direction, and omitting terms exponentially suppressed in t: 
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(4.16) 



* Here we omit the correct relativistic normalization of d since, as we will discuss later, 
we are not able to determine the polarizability in regimes where relativistic effects become 
sizable. 



where Ep = \Jp^ -\-M'^ — p'h- A comparison with Eq. (4.12) (non-relativistic approxi- 
mation and neglect of the motion of the hadron from ti to ^2) shows several new effects. 
First there are some extra relativistic normalization factors such as the factor multiplying 
B? which corresponds to Lorentz contraction. Secondly, the term linear in t arises from 
the motion of the hadron and reflects retardation effects. The latter are important for 
the Compton-polarizability of a charged point-like particle [15,16]. Furthermore, there are 
relativistic recoil terms of order /E"^. Most importantly, however, the momentum p of 
the initial state no longer factorizes, which makes it cumbersome to extract rms radii, po- 
larizabilities, and form factors from lattice measurements. For example, the polarizability 
becomes momentum dependent. The reason is that aud is defined as the second derivative 
of the energy M of a particle at rest in a background electric field, in our case d"^ M / dSud^d, 
where £u and £d act on up and down quarks, respectively. Since the electric field is not a 
Lorentz scalar, hadrons with different momenta p experience different polarizabilities. The 
polarizability even depends on the angle between the momentum p and the electric field 
E (this dependence is contained in the p-dependence of dhn which we have suppressed) . 

Due to these problems we will not discuss further how to extract the polarizabilities 
in the general case. Rather, it is clearly preferable to use momentum projected sources — 
which is a possible extension of this work in some cases — to extract the polarizability, as 
in the static case, by integrating the exponential term. In practice this is again achieved by 
integrating the difference between Eq. (4.16) and the full expression for the time dependent 
density-density correlation function, Eq. (4.3). 

We have determined the momentum distribution C| of our bag sources numerically 
by calculating the overlap between a localized source and a momentum projected source. 
Typical calculations show m~'^{p^) = rn~'^^^C3^]f / ^^^C^^ to be of the order of ten 
percent or less [17]. These are small enough that using the measured one could subtract 
the terms linear and quadratic in t and extrapolate to t ^ 00 in order to eliminate the 
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exponentially suppressed dipole-dipole terms to obtain the isoscalar rms radii 



\ {{ri) + m 




(4.17) 



(^charge) +! (^charge )^ nucleon 



Note that we have restricted our attention to u-d correlations for the usual computational 
reason that all observables can be evaluated using propagators calculated from the bag 
sources. In the future, u-u and d-d correlations need to be considered as well to disentangle 
the electromagnetic observables in the proton and neutron [4] . 

V. LATTICE RESULTS 

Monte Carlo calculation 

Density-density correlation functions were calculated on 12^ x 16 and 16^ lattices using 
20 quenched SU{?>) configurations generated by the Cabillo-Marinari heat bath method 
[18] with coupling p = = 5.7. For convenience, we have used the value of the lattice 
spacing determined from the sting tension, a ~ 0.2/m with ~ 1 GeV, with the result 
that the spatial length of our largest lattice is 3.2 fm. Clearly, for the purpose of comparison 
with the tree diagram path integral, the precise value of the lattice scale is inessential. 

Propagators were calculated at the three values of n shown in Table 1, for which 
the lightest pion mass in approximately 340 MeV. All the results shown in this paper 
used distributed bag model sources [13,19] at the first and last time slices, with hard- 
wall boundary conditions for the fermions to prevent quark propagation across the time 
boundary. 

The gauge fields were fixed to Coulomb gauge on the source time slices, and the bag 
radius was set to Ifm. As shown in Ref. [13], the bag sources project onto the hadron 
ground state extremely effectively, providing in general clean signals and broad plateaus 
in the number of time slices away from the source. Density operators were averaged over 
the central four time slices. In order to calculate the effect of momentum projection as 



mentioned in Section III, calculations for mesons were also carried out using uniform wall 
sources instead of bag sources on one time boundary. This, of course, produced significantly 
less of a plateau region on the wall source side. Analogous calculations for the nucleon 
were too noisy to be useful because integration over the wall of a third propagator which 
is not tied to a density operator requires a higher degree of phase cancellation than our 
present statistics could provide. 

Image Correction 

As described in section II, the physical correlation function for an isolated hadron 
is determined by subtraction of the tails of the densities of all first images. The most 
desirable way to perform this subtraction would be to use the known asymptotic behavior 
of the correlator, and the tree-level path integral described in this work provides one very 
attractive means to obtain this behavior. In analyzing the lattice results, however, we have 
used a self-consistent phenomenological analysis which utilizes only the lattice data itself. 

We iteratively improve the image corrections by approximately correcting for images 
using an appropriately defined parametric curve, least-squares fitting the parameters of the 
curve to the corrected data, and iterating to self-consistency. In practice, this procedure 
always yields a smooth, universal curve at large distances, as will be shown in examples 
below, and is insensitive to the precise form of parameterization. Note also that the 
iterative procedure is robust and achievement of a universal curve is a strong consistency 
check of the calculation and parameterization. If the parametric curve were too low, the 
image correction would be too small, and the corrected data would then be too high and 
inconsistent with the curve. The parametric curve is thus always driven in the proper 
direction, and consistency between the final curve and all the corrected data requires 
correct parameterization of the asymptotic density over a wide range of r. The actual 
form of the asymptotic density we have used is the following [19] 

p{r) = pQe-"^i^-5(^)('"-^)("^2-™i) (5.1a) 
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where 

g{r) = {1 + e-^^"--^^)-^ (5.16) 

This curve smoothly joins exponential decay with mass mi for r <<C i? with exponential 
decay with mass m2 for r :» i? with the transition occurring in the vicinity of R over a 
range b~^. This form has the advantage that one can easily select physical starting values 
and has sufficient flexibility for all the cases of interest in this work. 

Results 

Figure 8 shows how the self-consistent image correction procedure works in the case 
of the pion density-density correlation functions. The uncorrected data for a 16"^ spatial 
lattice are shown in (b). By symmetry arguments, the highest points at r/a = 8, 8y/2 
and 8^/3 are high by factors of 2, 4, and 8 respectively, and the other data in this regime 
display comparable finite volume effects. The self-consistent fit to the image corrected data 
is shown in (c), and one clearly sees that all the data now lie on a single universal curve 
which is very well fit by the self consistent solid curve specified by Eq. 5.1. Note that this 
curve is strongly constrained all the way out to the corner of the unit cell corresponding 
to 8\/3 lattice units or ~ 2.8/m. This is the curve that was shown in Fig. 5 and agreed 
in detail with the tree diagram path integral result. For comparison, the uncorrected 
lattice data for a 12^ spatial lattice are shown in (a), and the self-consistent fit to the 
image corrected data is shown by the dashed line in (c). This lattice is sufficiently small 
that in addition to the sum of tails of first images, one also observes errors beyond Ifm 
from interactions with the images. Given the accuracy of the tree level path integral in the 
region of 1-1.5 fm, even this result on an unphysically small 2.4 fm lattice can be combined 
with the tree level path integral to accurately describe the whole correlation function. 

The results for the pion density-density correlation function at all three values of n 
given in Table I are shown in Fig. 9. The solid curves are the self-consistent fits to image 
corrected data discussed above and have previously been shown in Fig. 5 where they 
agreed well with the tree level path integral results. For clarity in this and subsequent 
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graphs, nearby image-corrected lattice data are grouped into bins, and data within each 
bin are combined to a single value by means of a statistically weighted average of both 
the ordinates and abscissas. The top plot for K2 is thus the binned version of plot (c) of 
Fig. 8. The main conclusion from the results of this graph is that the self-consistent fits to 
the image corrected data are accurately determined for all three values of the quark mass. 
In comparing the lattice results with the slopes for the tt and p masses, we note that as 
already seen from the tree diagram results, in this regime the slope has not yet reached 
the pion slope. 

Analogous results for the rho density-density correlation function are shown in Fig. 10. 
In this case, we note that although one can still obtain a meaningful self-consistent fit, the 
statistical errors in the lattice data are significantly larger than for the tt. The most striking 
result is that as the quark mass decreases, the slope approaches the slope governed by the 
pion mass. We already noted in the discussion of the tree diagram analysis that there 
is a TTUjp effective coupling which produces a pion exchange diagram, and subject to the 
limitations of the statistics, these results indicate that the coupling is so strong that it 
dominates the decay at the separations addressed in this work. 

Finally, the results for the nucleon density-density correlation function are shown in 
Fig. 11. In contrast to the meson results, the plateaus in the number of time slices away 
from the wall are not always as clear as we would like, and at present, we have no theoretical 
argument why this problem arises for the nucleon. For K2, the plateaus become ambiguous 
for r/a> 12 and, indeed, one notes that the last three points appear systematically high. 
The problem is so pronounced for K5 that we were not sufficiently confident of the result 
to present it. 

The shape of the nucleon correlation function is somewhat different from that of the 
mesons. Whereas for mesons, a roughly constant asymptotic slope sets in around Ifm, 
we note that only in the vicinity of 2 fm does the nucleon slope begin to turn over and 
approach the rho slope. Physically, it is plausible that the three-quark core of the nucleon 
is more extended than the quark-antiquark component of the mesons, so that when the 
finite size of the relevant vertices is taken into effect, the tree diagrams only begin to 
describe the asymptotic slope at substantially larger distances. 



VI. SUMMARY, CONCLUSIONS, AND OUTLOOK 

In this work we have analyzed and performed measurements of correlation functions to 
study the pion, the rho and the nucleon. Our most important results are the development 
and successful testing of a scheme to correct systematically for images on periodic lattices, 
and the development of a tree level path integral to determine the asymptotic behavior 
of correlation functions. Furthermore, we have clarified the relation between the second 
moment of density-density correlators and the rms radii of hadrons. 

We have shown that as long as the interaction between hadrons and their images 
in adjacent unit cells is not too large, the difference between lattice gauge theory results 
for the density-density correlation function on a periodic lattice and an infinitely large 
lattice is simply due to superposition of the tails of the same correlation function arising 
from mirror images. In a similar situation in a non-gauge theory, non-diagonal interference 
terms between wave functions in adjacent cells would also contribute to the density-density 
correlations. In a gauge theory such interference terms are proportional to Wilson loops 
enclosing the periodic lattice and thus have the topology of Polyakov lines which are 
strongly suppressed. We have thus developed a practical method to use the long distance 
behavior of the density-density correlations on an infinite lattice to correct for the images 
on a finite volume lattice, and have demonstrated its effectiveness in an explicit lattice 
calculation. 

QCD at large distances can be most effectively described using hadronic degrees of 
freedom. As a model for the long distance behavior we have employed tree level diagrams, 
including only the lightest hadron in each channel. For the vector current correlations 
measured in this work, this is equivalent to vector meson dominance. We have used this 
model to calculate the asymptotic behavior of the same correlators we have treated on 
the lattice, and observe good agreement between self-consistent fits to the lattice data and 
predictions based on our tree level hadron model. This model also makes it clear why the 
large distance behavior of density-density correlations is not a simple exponential: since 
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the hadron propagates between the two vector meson vertices, in addition to the mass 
scale corresponding to the vector meson mass, the mass scale of the hadron itself is also 
relevant. We have also shown that when the density-density correlation is measured at 
spatial separations that are comparable to the temporal extent of the lattice, there are 
further nontrivial r-dependencies arising from incomplete momentum projection due to 
the finite extension of the lattice in the Euclidean time direction. An important practical 
result of this work is that because the most important finite size effects in the region r > 1.5 
fm are now understood, one can be confident of future measurements of correlations on 
the present physical volumes. 

Nonrelativistic arguments suggest a direct relation between the second moment of the 
density-density correlators and the rms radius, as defined by the slope of the form factor. 
This result arises because, nonrelativistically, the center of mass of the quarks separates 
and thus the contribution to the density-density correlation from dipole transitions in the 
intermediate state between the two current insertions is itself proportional to the rms 
radius. Unfortunately, this argument does not apply when the quarks are relativistic and 
the relation between the second moment of the density-density correlation and the rms- 
radii contains an a priori unknown dipole-dipole term. 

In this work we have only considered density-density correlations measured at equal 
time. In future work, it would be desirable to perform measurements of density-density 
correlations at unequal times for several reasons. By measuring the dependence of the 
second moment on the time difference, one could extract the polarizability of the hadron 
which is physically interesting and is being studied experimentally [15,16]. In addition, for 
sufficiently large times between the two current insertions, the dipole-dipole contributions 
arising from excited intermediate states are suppressed and only the contribution from the 
rms radius survives in the second moment. 
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APPENDIX A 



Density-Density Correlations for 
General Currents and Dirac Particles 



When we derived the large t-behavior of the second moments of the density-density 
correlation functions we factored the current matrix elements into the invariant form factor 
and a kinematic piece {q — p' — p) 



(A.l) 



Of course, for other currents one must use different kinematic factors. For example, 

{p'\uiu\p) = Fu{q'')ip + p) (A.2) 



or 



{p'\uu\p) = F^{q') , 



(A.3) 



where and are the vector and scalar up quark form factors of the hadron. Using these 
currents, the general structure of the terms appearing in the moments of the correlation 
functions will be very similar to Eq. (4.16), although the coefficients will be different. For 
example, for t = {t2 —ti) 
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where as usual, terms exponentially suppressed in t have been omitted. Since the Cps, 
as well as the isoscalar rms Eq. (4.17), can be determined independently, one can actually 
test Eqs. (A. 5) and (A. 6) by measuring Qa and independently. This provides a test of 
the Lorentz invariance which we have used throughout the derivations. Due to the coarse- 
grained structure of the lattice at short distances, Lorentz invariance is not guaranteed and 
it will indeed break down at some point for sufficiently large momenta. Testing Eqs. (A. 5) 
and (A. 6) is therefore useful to justify the Lorentz invariance assumption made in the text 
for "typical" momenta. 

We now consider the case of nucleon form factors. Here we will assume p = in the 
initial state. In practice this means that one neglects terms of order p^/M^ in the form 
factor. Since typical nucleon momenta in a lattice calculation are non-relativistic, this is 
a reasonable approximation. In fact, for non-relativistic nucleons p^/2M is of the order 
1/T. Hence //M^ is of the order 2/MT 0.1, where w e used M 1.2 a-^,T = Wa as 
typical values. 

Introducing invariant form factors via Refs. [20], [21] 



ip' \u^^u\ p) = u{p') 

and using 



u{p) (A.7) 



F-{q^) = Fi"(0) + i?.'^ + O^q") (A.8a) 
^^u = ^[FnQ) + F^m , (A.8b) 
as well as similar expressions for down quarks, one obtains for spin- 1/2 particles 

V 

-2^e-(^"-^'^)*dL4 • 

n 

(A.9) 
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Note that for 7° currents the magnetic matrix elements in Eq. (A. 7) are aheady of order 
q , i.e. they do not contribute to inelastic transitions. Hence the polarization terms in 
(A. 9) are purely electric and their time integral yields the u-d matrix element of the electric 
polarizability. 

The magnetic moments enter the expression for the second moment similar to the 
ground state rms radii. This is reminiscent of the Sachs form factor where also linear 
combinations of Fi and F2 appear [21]. The term linear in t enters with the same coefficient 
as for spinless hadrons. This should not be surprising, since it arises from the classical 
acceleration of the charged particle in an external electric field. Similar to the mesonic case, 
this term even survives for point- like particles. For Dirac particles (i?^ = 0, = 1/2M) 
the result still differs from point-like scalars and a finite "effective rms" — arising from 
the spin — remains. 

Of course, the non-relativistic approximation in (A. 9) is not really necessary, although, 
as we discussed in the context of spinless fields, the polarizability becomes rather difficult 
to determine. For the second moment, again omitting terms exponentially suppressed in 
t, one finds 



The only new term in (A. 10) is the term quadratic in the magnetic moment — indicating 
a double spin-flip contribution to the correlation function. 





(A.IO) 
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FIGURE CAPTIONS 

Fig. 1: Density-density correlation for non-relativistic positronium in QEDi+i on a lattice 
with N = 10 sites and periodic boundary conditions. The heavy solid curve shows the 
periodic solution on a circle and the dashed curve denotes the image corrected result. 
For comparison the result for an open line {N oo) is shown by the light solid curve, 
and this solution displaced by ±L is shown by the dotted curves to indicate the degree 
of overlap. 

Fig. 2: Fourier transforms of the density-density correlation functions shown in Fig. 1. The 
solid points denote the Fourier transform of the periodic density at the discrete lattice 
momenta = and the dashed curves show the Fourier transform of the solution 
on the open line. The solid curve in (a) shows the result of using the asymptotic form 
Eq (2.18) with b/a chosen to make the second derivative continuous, and the solid 
curve in (b) shows the errors introduced by a discontinuous second derivative arising 
from Eq. (2.16). The solid curves in (c) are obtained as in case (a) from an ensemble 
of ten data sets in which the periodic solution has been modulated with 1% relative 
random noise. 

Fig. 3: Time histories contributing to meson density-density correlation functions with peri- 
odic boundary conditions in the spatial direction. The sources creating the meson are 
denoted by the solid dots on the Euclidean time boundaries and the density operators 
are represented by the interior solid dots. Quark propagators on the front and back 
surfaces of the cylinder are denoted by solid and dashed wiggly lines respectively. 

Fig. 4: Space-time diagrams for {puPd)- The left sketch shows a typical time-history in the 
quenched approximation, with quark propagators denoted by wiggly lines. The sketch 
on the right shows the corresponding tree-level graph at large separations in which 
the pairs of qq propagators and the associated gluons are replaced by ground state 
meson propagators. 
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Fig. 5: Comparison of approximations to the pion density-density correlation function p(r) 
at large separation r. The sum of all tree- level diagrams without image corrections is 
given by the open squares. Image corrected results at the symmetry points are shown 
by solid octagons and interior points with negligible image corrections are indicated 
by solid squares. The best fit to the image corrected lattice QCD results is denoted 
by the solid curves and for comparison, the pion and rho slopes are shown by dotted 
lines. The dashed lines indicate the corrected asymptotic behavior obtained by taking 
the T — i> oo limit of the tree diagrams as discussed in the text. The three values of 
the hopping parameter k correspond to the masses given in Table I. The distance is 
calibrated both in lattice units (lower scale) and fm assuming a=0.2fm (upper scale). 

Fig. 6: Space-time diagrams for the pion (puPd) in the stationary tree-level approximation. 
The upper and lower graphs correspond to small and large spatial separations respec- 
tively. The left graphs correspond to the configurations in which the summand in 
Eq. (3.4) is maximal and the right graphs represent the non-summed geometry. Note 
the appearance of a p propagator in d. 

Fig. 7: Comparison between the stationary tree-level approximation to the pion (puPd) in the 
unprojected (full line) and the projected (dashed line) cases for relative separation r. 
For reference, the slopes corresponding to ttt-tt = 0.33, and nip = 0.615 are given by 
the dotted lines. 

Fig. 8: Lattice results for pion density-density correlation functions with and without image 
corrections. Uncorrected results with statistical errors on 12^ and 16^ spatial lattices 
are shown in (a) and (b) respectively, shifted by a decade for clarity. Shifted by an 
additional decade, the self-consistent fit to the asymptotic density for the 16^ lattice 
using Eq. (5.1) is shown by the solid line in (c), along with the image corrected lattice 
data using this fit. For comparison, the dashed line in (c) shows the analogous self- 
consistent fit obtained from the (12)^ data. All results are for the hopping parameter 
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K2 of table I. The radial distance is shown in lattice units on the bottom scale and in 
fm on the top scale. 

Fig. 9: Image corrected pion density correlation functions at the three values of k given in 
Table I. For clarity, the data have been aggregated into bins, and each correlation 
function is shifted by 3 decades. As in Fig. 8, the solid curves denote the self-consistent 
fits and the error bars denote the image-corrected lattice measurements. For reference, 
the slopes corresponding to m-j^ and nip given in Table I for the corresponding n are 
denoted by dotted lines. The radial distance is shown in lattice units on the bottom 
scale and in fm on the top scale. 

Fig. 10: Image corrected rho density-density correlation functions. The presentation and no- 
tation are the same as in Fig. 9. 

Fig. 11: Image corrected nucleon density-density correlation functions. The presentation and 
notation are the same as in Fig. 9. 
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TABLES 

Table I. Hadron masses at the three values of the hopping parameter k 
used in the lattice calculations. For reference, a bare quark mass = ^ ~ 2^ 
is also tabulated. The inverse lattice spacing is = 1 GeV. 







ruq 

(MeV) 


(MeV) 


rup 

(MeV) 


(MeV) 




0.16 


175 


691 (3) 


813 (4) 


1321 (10) 




0.1639 


95 


511 (5) 


698 (5) 


1097 (11) 


K5 


0.167 


40 


340 (7) 


615 (6) 


915 (16) 
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